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When a suppressed-carrier signal is squared, the carrier reappears in doubled form. An 
open-loop receiver can be used to deliver a recording of a band-limited waveform 
containing this carrier, whose amplitude and phase can be tracked by the radio science 
experimenter. 


I. Introduction 

In conventional radio science, a spacecraft transmits a 
sine-wave signal through a medium, which perturbs the signal 
in phase and amplitude. Inferences about the medium are 
made from observations of these phase and amplitude varia- 
tions. Because the present DSN telemetry system uses a 
residual carrier, open -loop or closed-loop radio science can be 
carried out in the presence of telemetry modulation. Figure 1 
shows an open-loop processor. 

The multimegasymbol telemetry system now under devel- 
opment can use a suppressed-carrier signal having no discrete 
spectral component at the carrier frequency. Doppler phase 
and signal amplitude can be extracted by a Costas-loop 
receiver; hence, closed-loop radio science can be carried out as 
before. In fact, as Ref. 1 shows, this scheme enjoys a 15- to 
30-dB advantage over the residual-carrier scheme (with an 
80-deg modulation index and no subcarrier) because the latter 
uses only 1/33 of the total signal power. 

Our purpose is to ask whether open-loop radio science can 
also be carried out with a suppressed-carrier signal. After some 
tentative investigation into the structure of the optimal 
receiver for extracting the phase and amplitude of the 
suppressed carrier, we decided to take a more practical turn. 


Since the Costas loop has the same performance as a squaring 
loop (Ref. 2), why not use an open-loop heterodyne receiver 
to process the square of the residual-carrier waveform? Such a 
processor is shown in Fig. 2. The recorded waveform contains 
a doubled version of the carrier. 

The comparison between this scheme, which we call 
squared-suppressed-carrier (SSC), and open-loop reception of a 
residual carrier (RC), almost duplicates the closed-loop com- 
parison in Ref. 1 . The reason for this is that the tape recording 
of the squared waveform is analyzed by a phase-locked loop in 
computer software to extract the phase and amplitude of the 
recovered doubled carrier. In effect, this is nothing but a 
squaring loop. Nevertheless, we have taken some care to 
examine exactly what information (and noise) is contained in 
the square of the suppressed-carrier waveform, as influenced 
by the dispersion of the medium and of the presquaring filter. 

The comparison between RC and SSC is based on the noise 
variances of the measurements of the phase and amplitude of 
the perturbed carrier. In the SSC case, there is, of course, no 
discrete carrier for the medium to perturb, yet we shall see 
that the reconstructed doubled carrier in the squared wave- 
form contains almost the same information about the medium 
as the perturbed residual carrier. In fact, if the dispersion of 
the medium over the bandwidth of the suppressed-carrier 
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signal can be neglected, then the S SC receiver of Fig. 2 delivers 
a signal whose phase is twice the RC phase and whose 
amplitude is the square of the RC signal amplitude. Therefore, 
the experimenter must halve the phase and take the square 
root of the amplitude to obtain phase and amplitude informa- 
tion equivalent to what the RC setup gives. The effect of 
dispersion is assessed in Section V. 

We shall show that the SSC scheme enjoys a considerable 
signal-to-noise-ratio (SNR) advantage over the RC scheme for 
the measurement of phase and amplitude. The results can be 
summarized here. The RC and SSC output tapes (Figs. 1 
and 2) are both analyzed by phase-locked loops in software to 
recover phase and amplitude; let the loop bandwidth in both 
cases be B L . Assume that both the RC and SSC received 
waveforms (the inputs in Figs. 1 and 2) have total signal power 
A* and noise spectral density Vi N 0 . Denoting the estimated 
carrier phase by \p e , we have 


A. Transmitted Signals 

These are 


y/2 COS 0 COS CO / 

T C 

(residual carrier) 

0) 

\fl d(t) cos co.r 

(suppressed carrier) 

(2) 


where 6 = modulation index, d(t) = data modulation, modeled 
by independent random ±l’s with a symbol time T. In our 
situation, 0.1 MHz < 1/T < 30 MHz. The spectral density of 
d(t) is 

S d (K) = Tsinc 2 04X7) (3) 

The random Fourier expansion of the stationary process d(t) 
can be written symbolically* as 


1 

o 2 ($ e ) 


N o B l 


cos 2 e 


(residual carrier) 


d(f) = 



D(X)e iX, d\ 


(4) ‘ 


1 € 

= /?' (squared suppressed carrier) 

o 2 ($ e ) N o B l 

For RC, 0 is the modulation index. For SSC, 0' is a 
degradation factor, given in Eq. (30) below, caused by the 
distortion of the presquaring filter G ( fcj) (Fig. 2) and the noise 
induced by the squaring operation. These equations are also 
found in Ref. 1 , which compares RC and Costas-loop tracking 
in the multimegasymbol telemetry environment. For 6 = 80 
deg and a symbol SNR between 0 dB and 10 dB, we have 

cos 2 0 = -15.2 dB, -3.1 dB<0' <-0.8 dB 

for a proper choice of the bandwidth of the presquaring filter 
(see Fig. 3). It follows that SSC has an advantage of 12 dB to 
14 dB over RC for the measurement of phase. The same 
conclusion holds for the measurement of log-amplitude. This 
advantage can be used either to decrease the total noise 
variance of the measurement or to increase the" measurement 
bandwidth B L . 

II. Models 

We shall set up conceptual models for open-loop reception 
of RC and SSC signals. To be most fair to the RC setup, we 
shall assume that a subcarrier is used. In what follows, to and X 
denote radian frequency. 


B. External Medium and Received Waveforms 

This is what the radio science experimenters want to probe. 
We represent it by a time-varying transfer function G e (u, t ), 
the “external filter.” It operates as follows: Any RF exponen- 
tial Ae iu>t emerges from the medium as the narrowband 
waveform G e ( co, t)Ae iu> r , and the principle of superposition 
holds for linear combinations of complex exponentials. Fur- 
ther, G e ( co, f) is scaled so that it accounts for the power 
transmitted by the spacecraft, the space loss, and the gains of 
the transmitting and receiving antennas. Thus, the received 
residual-carrier waveform is 


~” 7 F”COs 6 G e ( co,, t) e“ c ‘ + c.c. + njt) 


= vTcos e A e {t) cos (to c t + 4> e (t)) + n w (t) (5) 

where c.c. means complex conjugate, G e (u> c , t) has amplitude 
A e (t) > 0 and phase i// e (r), and « w (f) is white Gaussian noise 
with spectral density Vi N 0 . We presume that the desired 
radio science information is contained in ^4 e (0 and i^ e (r). 


* Actually, D(\) d\ represents a process with orthogonal increments. In 
the Appendix, we shall approximate d(t) by a process with a discrete 
Fourier expansion. 
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Applying G e to Eq. (2), with d(t) given by Eq. (4), we 
obtain the received suppressed-carrier waveform as 


The output of the internal filter in Fig. 2 is jc(/) =s(t) + n(t), 
where 


— —d ( t)t ,u>cr + c.c. 

y/2 e 


( 6 ) 


where 


■I 


iXf 


d e (t) = I G e (<o c + X,r)D(X)e‘ ‘ d\ 


(7) 


s(r) = — —d r (t)e ,u>cT + c.c. 

V2 


n(t) = gft)*n w (t) 

and g.(t) is the impulse response of G.. Then 

x 2 (t) = s 2 (r) + 2s(t)n(t) + n 2 (f) 


(9) 

( 10 ) 

(11) 


is the data waveform after distortion by the baseband filter 
G e ( co c + X, t ). If 

G e^c ' X - 0 = G~ e ( w , ’c + X - f ) 

were to hold (the bar means complex conjugate), then d e (t) 
would be real. In general, d e (t ) is complex. 

C. Open-Loop Receivers 

Stripped to its conceptual essentials, an open-loop receiver 
for the residual-carrier waveform Eq. (5) is shown in Fig. 1. 
The heterodyne frequency f h is chosen such that 0 < f c - f h 
<B V , where B v is the “video” bandwidth, ranging from 100 
Hz to 50 kHz in the DSN Multi-Mission Receiver (Ref. 3). The 
video waveform is sampled at 2 B v and recorded. 

To process the suppressed-carrier waveform Eq. (6), the 
receiver of Fig. 2 is suggested. Before being squared, the 
waveform passes through an “internal” bandpass filter G ( - with 
a nominal passband f c ±B. We shall show that the squared 
waveform contains a discrete component at 2 f c containing the 
radio science information. Mixing this with 2 f h gives another 
video waveform with bandwidth 2 B v . This processor will be 
called the squared-suppressed-carrier receiver. 


111. Analysis of the Squared-Suppressed- 
Carrier Receiver 

Let G( co, t) = G g ( co, t)G.(u>), the overall filter. In analogy 
with Eq. (7), set 


d G (t) = \ G(u> c + X, t)D(K)e iXt d\ (8) 


Let us first examine the S X S signal s 7 (t). In the Appendix 
it is shown that s 2 (t) is the sum of two components, one with 
discrete spectral lines and the other with a spectral density. 
The discrete component has spectral lines at frequencies 2 f + 
n/T ( n = 0, ±1 , ±2, • • ■) plus a DC term that we ignore. 

The central component at 2 f c , the reconstructed doubled 
carrier, is the narrowband waveform 

1/2 F(0, t)e‘ 2uict + c.c. (12) 

where 

/* oo 

no, 0 = G(w c + X, t)G(u e - X, t) sine 2 (14 XT) ^ 

*'-00 

(13) 

(a special case of the function Y (co, t) defined in the 
Appendix). If G(w, t ) were equal to G e (u> c , t ) for all u>, then, 
because sine 2 Q/AT) integrates to 1, Eq. (13) would become 
T(0, r) = G 2 (to c , t). This means that Eq. (12) would carry the 
same radio science information as the residual-carrier signal in 
Eq. (5) except that here the carrier would appear in doubled 
form; Eq. (12) would be 

A]{t) cos (2co c t + 2i p e (t)) (14) 

In practice, G will have dispersion due to the external medium 
and the receiver, especially the ±B bandpass. The actual Eq. 
(13) gives T(0, t) as a weighted average of G(co c + X, t)G((x> c - 
X, /). We shall say more about this in Section V. 

The spurious S X S components at 2 \f + n/T {n 4= 0) will be 
ignored because we assume that the video bandwidth 2 B v is 
less than 1/7'. As with Costas loop tracking, one must be 
careful not to mistake a spur for the main component (Ref. 4). 
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The continuous-spectrum component of s 2 (t), the S X S 
noise, is almost small enough to be ignored in this study. 
Nevertheless, we shall set down its spectral density here for the 
record. Let co, = 2-rr/T. Define the function 


S (2co , t) = 'AN 2 

nn y c’ ' 0 



I G.(w c + X) 


X G,(« e - X)| 2 


dX 

2tt 


( 19 ) 


4>(0,X, t) = ^ G(co c + X + ko> x , t)G(oJ c - X - fcojj, t) 

k--°° 


(Because n(t) is Gaussian, the spectral density of n 2 {t) is 
2S„(o>)*S n (<o).) 


X sine 2 (A (X + k(o i )T) (15) 

(a special case of the function 4>(to, X, /) defined in the 
Appendix). This function (of X) is u>, -periodic, and its average 
over a period is just y(0, t). The spectral density of the S X S 
noise at 2 f c is now 


1^(0, X, f)- y(0, r)| 2 ™ 

u 

(16) 

which is 772 times the mean-square variation of <1> over a 
period. This is a time-dependent, short-term spectral density, 
valid for integration times that are short compared to the 
fluctuations of G e ( to, r). 


IV. Comparison of Open-Loop Receivers 

A. Ratio of Power to Noise Spectral Density 

Before proceeding to the comparison of the noise variances 
of phase and log-amplitude, it is convenient to compute 
power/noise ratios for the RC and SSC receivers. For the RC 
receiver, the ratio of signal power to one-sided noise spectral 
density (NSD) of the waveform Eq. (5) is 

A 2 

=~ cos 2 6 (residual carrier) (20) 

We have dropped the dependence of A e on time. 

To get a simple expression for the squared-suppressed- 
carrier (SSC) receiver, let us assume that 


Equation (16) is equivalent to Eq. (17) of Ref. 5. 

G e (w c + X, r) = G e (w e , f), (IXI < 2nB) (21) 

The spectral densities of s(r) and n(t ) in Eqs. (9) and (10) 
for frequencies near f c are C ; (w c + X) = 1, (|X| < 2nB) 

( 22 ) 

S/cu. + X, r) = 561 G(w c + X, t)\ 2 T sine 2 (A XT) = 0> otherwise 


= KJVolG^+l)! 2 (17) 

Therefore, the spectral densities at 2f c of the S X N noise 
2 s{t)n(t) and the NX N noise n 2 {t) are (Ref. 6) 


This means that we are neglecting the dispersion of the 
external filter over the passband of the internal filter, which is 
assumed ideal. Then the doubled reconstructed carrier, Eq. 
(12), is 


5 «( 2w c- 0 



I G(w c + X, r) G£fti c - X)| 2 


A] (/) a cos (2 coj + 2\p e (t)) 


where 


(23) 


X sine 2 (A XT) 


TdX 

2 it 


(18) 


a = 


I 


2 nBT 


-2rr BT 


sine 2 '(A x) 

2.71 


(24) 
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This a (the same as in Ref. 1) measures the signal loss due to 
the distortion of d(t ) by the receiver. 

The noise spectral densities from Eqs. (18), (19), and (16) 
become (with time dependence dropped) 


For SSC we must remember that the signal being tracked is 
the doubled carrier Eq. (23) (heterodyned to video). Thus, the 
usual linear-loop phase error variance is really the variance of 
2\p e . Computing the variance of i p e introduces a factor of four 
that cancels the 1 /4 in Eq. (26) and gives 


N 0 A > 

(SX N) 

1 

A l , 

n 2 0 b 

(NX N) 

(25) ° 2 ti e ) 

= 6 (squared-suppressed-carrier) (29) 

N o B l 

UA*Th 2 

(S X S) 

where, as in Ref. 1 (also see Refs. 2 and 5), 


where p 2 , a function of BT, is the mean-square oscillation of 
over a period 2n/T. In the cases given below, the S X S 
density is almost negligible. 

Now let R = A 2 T/N Q , the symbol SNR. Equations (23) and 
(25) give 


ft = 


Ra 2 

Roc + BT 


(30) 


The S X S noise, small to begin with, has been suppressed by 
the coherent phase detection. 


Power 

NSD 


HA* a 2 
e 


2(N 0 A 2 e a+N 2 n B'+KA*TnJ 


e ^2 J 


N n 


Va y (squared-suppressed-carrier) (26) 


where the degradation factor y is given by 


y = 


Ra 2 


Ra+BT+ HR 2 n , 


(27) 


which is the same as 6' of Ref. 1 (see Eq. (30) below) except 
for the small S X S noise term H R 2 U 2 . The next subsection 
will show why the factor 1 /4 in Eq. (26) is not counted as a 
degradation. 


Because the SSC receiver does not have to deliver telem- 
etry, the presquaring filter bandwidth 2 B can be chosen to 
minimize the phase noise variance. It can be assumed that the 
symbol rate 1 jT is set so that the nominal symbol SNR R = A 2 
T/N 0 has some reasonable value, say between 0 dB and 10 dB. 
Given R, one can choose BT to maximize ft. (Recall that a is a 
function of BT from Eq. (24).) The resulting BT and 6' are 
given as functions of R in Figs. 3 and 4 for -4dB</?< 
12 dB. (At a value of R slightly more than 12 dB, the optimal 
BT jumps discontinuously to a value greater than 1 .) Also 
shown is the residual carrier suppression cos 2 9 for several 
values of the modulation index 6; the ratio ft /cos 2 9 is the 
advantage of SSC over RC. 

Let the amplitude of the signal in the video waveform be 
extracted by a coherent amplitude detector (Fig. 5) whose 
low-pass filter has one-sided noise bandwidth B A . With perfect 
phase tracking, the noise variance of the amplitude estimate A 
is 


B. Amplitude and Phase Noise Variances 


o 2 (A)=N a B a (31) 


Consider now the actual quantities of interest, namely A e 
and \fr e . Let us suppose that for both RC and SSC the phase of 
the signal in the video waveform is tracked (in nonreal time) 
by a phase-locked loop with one-sided bandwidth B L . For RC 
(the waveform of Eq. (5) heterodyned to video), the variance 
of the loop phase estimate i p e is given by 


where N A is the one-sided spectral density of the baseband 
noise from the amplitude detector. 

For RC we have A = A e cos 9, N A =N 0 ; hence 



n o B l 


cos 2 9 


(residual carrier) 


(28) 



a 2 (A e ) 


N o B a 


cos 2 9 


(residual carrier) 


(32) 
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For SSC we have A = (a/\f2) A 2 , N A = 2 N Q A 2 a + 2A r 2 B + 
2AgTfx 2 - Since o(A)/A ~ 2a(A e )lA e , the % factor in Eq. 
(26) again goes away, and 



a 2 (A e ) N o B a Ra + BT + R 2 ii 2 


TdX 


X G .( co - X) sine 2 (% X T) ^ 

1 c Z7T 


4(2 h -h]) r°° 

+ — J^,.K + X) 


(squared suppressed carrier) (33) 

The degradation factor in Eq. (33) is the same as y (Eq. (27)) 
except that the small S X S noise term has been doubled by 
the coherent detection. It can still be ignored in this context. 
Again, one can use Fig. 4 to compare SSC to RC. 


V. The Effect of Dispersion 

The suppressed-carrier signal, although its bandwidth may 
be as much as 60 MHz, is still a narrowband signal because the 
carrier frequency is at least 2.3 GHz. Dispersive effects of the 
medium may therefore be negligible. Nevertheless, let us 
measure the dispersion of the external filter G e (cj, t) near co c 
by the function 


H(K t) = 


G e ( w c + x - * ) 

G e ( co c , t) 


(34) 


X G.(oj c - X) sin 2 (^ XT) ^ 


( 36 ) 


The second term carries the dispersive effect. Although one 
cannot measure the dispersion by observing Y (0, /), one could 
use prior bounds on X, and h 2 , and their variation with time, 
to estimate how much the ratio Y (0, t)jG 2 (c o c ,t) varies with 
time. 


VI. Conclusions 

By squaring a suppressed-carrier signal, it is possible to 
provide the radio-science experimenter with a tape containing 
the reconstructed doubled carrier. For measuring amplitude 
and phase, this scheme has typically a 12-dB to 14-dB 
advantage over the present residual-carrier scheme. 


Then G( gj c + X, t) = G e (w c , t ) G i (co c + X)H(\ t). If H (X, t) 
does not depend on t, then by Eq. (13) the complex amplitude 
Y (0, r) of the reconstructed doubled carrier is still propor- 
tional to G\ (w c , t). In any case, assume that for gj c + X in the 
passband of G i we have the development 

H(X,t) = l+h i (t)\ + h 2 (t)\ 2 + •■: (35) 


If the dispersion of the medium over the presquaring filter 
passband is neglected, then the reconstructed doubled carrier 
yields the same information as the residual carrier; one merely 
takes the square root of the amplitude and halves the phase. 
There is a phase ambiguity of n instead of 2rr. This is of no 
consequence if the phase is tracked continuously. 


Assume that G f (co, t) is known. Neglecting the remainder in 
Eq. (35), we have 


no, t ) 


G](u c , t ) 



G,(oj c +X) 


It is anticipated that the effect of media dispersion on the 
doubled carrier reconstruction will be small because the 
bandwidth of the suppressed-carrier signal is small compared 
to, for example, the difference between Sband and Xband. 
The effect of this dispersion can be estimated. 


91 



References 


1. Lesh, J. R., ‘Tracking loop and modulation format considerations for high rate 
telemetry,” DSN Progress Report 42-44, Jet Propulsion Laboratory, Pasadena, 
California, 1978. 

2. Lindsey, W. C., and Simon, M. K., Telecommunication Systems Engineering , Prentice 
Hall, 1973. 

3. Donnelly, H., and Nishimura, H., “Multi-Mission Receiver (MMR),” DSN Progress 
Report 42-52, Jet Propulsion Laboratory, Pasadena, California, 1979. 

4. Hedin, G. L., Holmes, J. K., Lindsey, W. C., and Woo, K. T., ‘Theory of False Lock in 
Costas Loops,” IEEE Trans. Commun., COM-26, 1-12, Jan. 1978. 

5. Simon, M. K. and Lindsey, W. C., “Optimum Performance of Suppressed Carrier 
Receivers with Costas Loop Tracking,” IEEE Trans. Commun., COM-25, 215-226, 
Feb. 1977. 

6. Davenport, W. B., Jr. and Root, W. L., An Introduction to the Theory of Random 
Signals and Noise, McGraw Hill, 1958. 



f 

c 


RESIDUAL-CARRIER 

WAVEFORM 



f 


TAPE 
TO USER 


Fig. 1. Conceptual model of residual-carrier, open-loop receiver 
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Fig. 3. Finding the optimal pre squaring filter bandwidth 2S for the 
squared-suppressed-carrier (SSC) receiver [T = symbol time) 
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Fig. 4. Comparative losses of the SSC and residual-carrier (RC) 
receivers for several RC modulation indices 6 
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Fig. 5. Coherent amplitude detector 
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Appendix 

S x S Components 


Our aim here is to compute the discrete -spectrum and 
continuous-spectrum components of s 2 (t), where s(t) is given 
by Eq. (9). To do this, we have found it expedient to 
approximate d(t) by a waveform d N (t) of period NT, defined 

by 


d N (t) =d. (jT<t<(j+l)T),0<j<N- 1) 


d N (t + NT) = d N (t) (ail t ) 


(A-l) 


where d 0 , • ■ • , d N _ , are independent random variables taking 
the values 1 and -1 with probability 1/2. Then d N (t) has the 
discrete Fourier expansion 


VO = Z 


(A-2) 


where the X spacing is AX = 2n/(NT). Results will be obtained 
by letting N~*°°. 

Let d(t) be replaced by d N (t). Ignoring DC components, we 
have 

; 2 ui „ r 


s 2 (f) = Vie c y(t) + c.c. — . 


where 


AO 


[? 


G { co + X, t)D(K)t 




= L nu.tW'** 


(A-3) 


co 


The co spacing is Ac o = 2vj(NT). 
From Eq. (A-3), 


To compute the mean and variance of Y (co, t ) we need the 
second and fourth moments of the D (X) process. Let co, = 
2rr/7’. From Eqs. (A-l) and (A-2), 


Z>(X) 


p -i\T , , W-i 

e - i i , -//xr 


-/XT N 




1=0 


(X = ku>JN,- «><*:<«>) 


Then 


^[^(Xj)/)^)] = -^ sine (^XjT) sine 


(X, + X 2 = fcco,) 

= 0, (X, + X 2 it 0 mod co,) (A-5) 

E[D(\ t )D(X 2 )D(\ 3 )D(X A )] 

= EiD^DfX^EiD^D^)] (parti) 

+ £-[Z)(X 1 )T»(X 3 )]£'[D(X 2 )r>(X 4 )] (part 2) (A-6) 

+ f[D(X 1 )i)(X 4 )]T[2)(X 2 )D(X 3 )] (part 3) 

+ part 4 

part 4 = (-l) k sine (& X T) • • • sine (& X T), 

/V 3 

(X, + X 2 + X 3 + X 4 =ku> x ) 

= 0, (X, +X 2 +X 3 +X 4 ¥=0modoo,) 

If d{t) were a Gaussian process, then part 4 would be zero. 
From Eq. (A-4), 


F(co, 0 = Z + ^0C(w c + w-X, f) 


co 


£T(co, r) = 2 G (co c + 0G(w c + co - X, f) 


X Z?(X)D(co - X) 


(A-4) 


F[Z>(X)D(co - X)] 
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£\Y(u,t)\ 2 = £ £ G(w c + X,OG(« e +w-X,0 
x x’ 

• G(co c + X', t)G(u) c + co - X', /) 

• £ [D (X)D (w - X)Z>(- X')D(X' - co)] A-7) 

where the bars mean complex conjugate. It is convenient to 
use the notation 

F(co, X, t) = G(co c + X, r)G(co c + co - X, /) 

X sine (Yi XT) sine (14 (co - X) F) 

(A-8) 

in terms of which the expectation of Y (co, t) is: 


EirCiKOj.OI 2 = bEY^co ,,?)! 2 

which says that the variances of the real and imaginary parts of 
Y(nu v t) are zero, and hence 

F(mo, , f) = (- 1)" J F(nu> r X, t) ^ (A-9) 


a nonrandom quantity. When n = 0 we get Eq. (13) for the 
complex amplitude of the recovered doubled carrier; when 
n#0we get the spurs, which are expected to be smaller than 
the main peak because the functions sine ( Yi (X- nco,) T) 
are orthonormal. These discrete components necessarily have 
nonrandom amplitudes, otherwise s 2 (t) would have strange 
nonergodic properties. 


EY(nu v t)-^~ £ ^("Wj.X.0 (« = 0, ±1 , • • •) 

x 

£Y(u,t)= 0 (co^Omodojj) 

When Eq. (A-6) is used to compute the fourth moments in Eq. 
(A-7), the variance of Y (co, t) takes the form 


Let co # o mod co, . Then EY (co, t) = 0 and ElY( co, t)| 2 
= 0(1 IN). This means that as N -*■ °° the discrete spectrum 
near co tends to a spectral density. We can let 1/N = 
T Aco/(27r). Then 


S (2co„ + co, t) =1/4 Urn 
Aur-*-0 


£|y(co, r)| 2 
Ao)/(2tt) 


4 

e i y(co, t) i 2 = £ P art 1 

1=1 


where 

part 1 = |£T(co,OI 2 


part 2 


part 3 — 
N 2 


E 


£ F(co, X, t)F(co, X + fccjj, r) 

fr = -~ 


Define 


T 

2 



F( co, X.0 £ F(co, X + fccOj , r) 

fe=“> 


rax 

2n 




£(co, X, t) 


rax 2 

2 ir 


(A-10) 


part 4 = - 


2 

W 3 



2 


*■(«, X, 0 


<l>(co, X, t) = £ F(co, X + fccOj , r) 

k =- oo 


We are going to interpret these by letting When we 

do this, all the sums become integrals in which l/N becomes which is an coj-periodic function of X. In terms of it, Eq. 

TdXlQ.it). Parts 2, 3, and 4 are all O(l/A0- For co = ncoj we (A-10) becomes 

have in the limit 


EY(ntJ v t) = (-1)" 
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. . TdX 


2 


\*(cj,\t)-Z(cj,t) | 2 ^ 
(A-l 1) 


where 


Z(co,t) = f F(u.\t) ~ = f ™ 

J -oo •'o 

(A- 12) 

Note that Z(nco,, t) = (- 1)" V(nu) 1 , t) (n= 0, ±1, • • •). The 
special case to = 0 gives Eq. (16) of the text for the S X S noise 
spectral density at 2 / 
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